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ABSTRACT 

Context. The inner disc, linking the thin disc with the bulge, has been somehow neglected in the past because of intrinsic difficulties in 
its study, due, e.g., to crowding and high extinction. Open clusters located in the inner disc are among the best tracers of its chemistry 
at different ages and distances. 

Aims. We analyse the chemical patterns of four open clusters located within 7 kpc of the Galactic Centre and of field stars to infer the 
properties of the inner disc with the Gaia-ESO survey idr2/3 data release. 

Methods. We derive the parameters of the newly observed cluster, Berkeley 81, finding an age of about 1 Gyr and a Galactocentric 
distance of ~5.4 kpc. We construct the chemical patterns of clusters and we compare them with those of field stars in the Solar 
neighbourhood and in the inner-disc samples. 

Results. Comparing the three populations we observe that inner-disc clusters and field stars are both, on average, enhanced in [0/Fe], 
[Mg/Fe] and [Si/Fe]. Using the idr2/3 results of M67, we estimate the non-local thermodynamic equilibrium (NLTE) effect on the 
abundances of Mg and Si in giant stars. After empirically correcting for NLTE effects, we note that NGC 6705 and Be 81 still have a 
high [tt/Ee]. 

Conclusions. The location of the four open clusters and of the field population reveals that the evolution of the metallicity [Fe/H] and 
of [a/Fe] can be explained within the framework of a simple chemical evolution model: both [Fe/H] and [a/Fe] of Trumpler 20 and 
of NGC 4815 are in agreement with expectations from a simple chemical evolution model. On the other hand, NGC 6705, and at a 
lower level Berkeley 81, have higher [a/Fe] than expected for their ages, location in the disc, and metallicity. These differences might 
originate from local enrichment processes as explained in the inhomogeneous evolution framework. 

Key words. Galaxy: abundances, open clusters and associations: general, open clusters and associations: Individual: Trumpler 20, 

NGC 4815, NGC 6705, Berkeley 81, Galaxy: disk 


1. Introduction 

Open clusters are among the best tracers of the chemical evo¬ 
lution of the Galactic thin disc (see, e.g., Friel 1995; Bragaglia 
& Tosi 2006) from the very distant outskirts (e.g. Sestito et al. 
2008; Yong et al. 2012) to the inner regions, close to the Galac¬ 
tic bulge (see, e.g., Magrini et al. 2010, 2014). In particular, the 
clusters in the so-called inner disc, i.e. the part of the Galactic 
thin disc located at Galactocentric distances <8 kpc', have been 
so far relatively little explored; indeed, few clusters are known 
and spectroscopically observed (e.g Carretta et al. 2005, 2007; 
Sestito et al. 2007; Magrini et al. 2010; Carraro et al. 2014a,b,c; 
Magrini et al. 2014). This is due both to observational limits, i.e., 
high reddening and crowding towards the Galactic Centre, as 
well as to high mortality of clusters in regions where the density 

* Based on observations collected with the FLAMES spectrograph at 
VLT/UT2 telescope (Paranal Observatory, ESO, Chile), for the Gaia- 
ESO Large Public Survey (188.B-3002). 

' In the whole paper the adopted Rgc.o is 8 kpc (see Malkin 2013) 


of stars is higher (Freeman 1970; van der Kruit 2002). However, 
this region is of great importance because it constitutes a con¬ 
nection between the properties of the bulge and of the thin/thick 
disc. 


The Gaia-ESO Survey (Gilmore et al. 2012; Randich 
et al. 2013) is providing high-resolution spectroscopy with 
FLAMES®VET and homogeneous stellar parameters and abun¬ 
dances for a large sample of Galactic open clusters. The sample 
includes clusters located throughout the disc, at different dis¬ 
tances from the Sun, and with ages ranging from their initial 
phases (with the exception of the embedded ones), i.e., a few 
Myr after their formation, up to several Gyr (Randich et al., in 
prep.). The selection of targets to be observed within clusters 
is as uniform as possible: briefly, Gaia-ESO observes unbiased 
samples of cluster candidates with GIRAFFE, i.e. their selection 
is based only on their location of the colour-magnitude diagrams 
and no other information, as radial velocity, is used to ensure 
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their membership, and more secure cluster members with UVES 
(see Bragaglia et al., in preparation for more details). 

Observations for several clusters were completed in the first 
18 months of the Gaia-ESO Survey. Analysis is now completed 
and the products have been released to the consortium as Inter¬ 
nal Data Release 2, hereafter idr2/3. The sample includes four 
very interesting intermediate-age/old clusters in the inner disc: 
NGC 6705, Trumpler 20, NGC 4815, and Berkeley 81. The first 
three have been discussed in a series of previous paper using 
Gaia-ESO idrI recommended values (e.g., Magrini et al. 2014; 
Tautvaisiene et al. 2015). In the present paper, we not only con¬ 
sider the updated idr2/3 analysis, with important differences in 
the determination of the final parameters with respect to idrI 
(see Smiljanic et al. 2014, for details), but we add to the sample 
Berkeley 81 (hereafter Be 81), which is the innermost cluster in 
the idr2/3 sample. We aim at comparing its abundance pattern 
with both field stars located in the Solar neighbourhood and to¬ 
ward the Galactic centre, and with the other three clusters located 
in the inner disc. 

Be 81 is an intermediate-age cluster located at approximately 
6 kpc from the Galactic centre (1=34.505, b=-2.068). It has 
been recently studied by Donati et al. (2014a), who, compar¬ 
ing observed and synthetic colour-magnitude diagrams, found 
an age close to 1 Gyr, a large reddening E(B - V) =0.9, a dis¬ 
tance modulus of 12.4, and a metallicity close to Solar. Pre¬ 
vious literature data on Be 81 are quite scarce: BVI photom¬ 
etry of a portion of Be 81 was obtained by Sagar & Griffiths 
(1998), from which they derived an age of ~1 Gyr and a redden- 
ing£’(B-y) =1. The only determination of the metallicity of this 
cluster is by Warren & Cole (2009) through calcium triplet spec¬ 
troscopy. They obtained a rather uncertain sub-solar metallicity 
([Ee/H]=-0.11+0.15 dex), somewhat unusual for its internal lo¬ 
cation in the Galactic disc. This is likely related to the limiting 
magnitude of Warren & Cole (2009) that did not allow them to 
reach stars as faint as the cluster giant branch and so they likely 
were not able to identify any of the cluster member stars. 

The paper is structured as follows: in Section 2 we describe 
the analysis of Be 81 in idr2/3. In Section 3 we derive its funda¬ 
mental parameters, whereas in Section 4 we compare the chemi¬ 
cal patterns of field and cluster populations. In Section 5 we dis¬ 
cuss the chemical patterns, and in particular the a-enhancement, 
in the context of our current knowledge of the disc’s evolution. 
In Section 6 we give our summary and conclusions. 

2. Sample and analysis of Be 81 

In Be 81, 14 stars were observed with UVES using the setup 
U580 (R=47000, AA 4800-5800 A); 266 stars were instead tar¬ 
geted with GIRAEEE using HR15N (R=17000, AA 6470-6790 
A) and HR09B (R=25900, AA 5143-5356 A) (16, 80, and 170 
stars observed with both setups, HR09B only, and HR15N only, 
respectively). In Eigure 1 (right panel) we show the colour- 
magnitude diagram of the inner 6 arcmin of the cluster. The 
UVES targets are candidate red clump giant stars, while the 
candidate main sequence stars were observed with GIRAEEE. 
In the upper left panel, the field of view of the photometric data 
in Donati et al. (2014a) obtained with the Large Binocular Cam¬ 
era at the Large Binocular Telescope (LBC@LBT) is presented. 
Einally in the lower left panel, we show the distribution of the 
radial velocities versus distance from the cluster centre. 

The spectrum analysis of the Gaia-ESO data is organised in 
several working groups (WGs) that provide astrophysical param¬ 


eters and element abundances. In particular, the analysis of the 
UVES and GIRAEEE spectra discussed in the present paper was 
performed by three WGs: one dedicated to EGK-type stars ob¬ 
served with UVES (WGll- see Smiljanic et al. (2014) for de¬ 
tails), one for the same type of stars targeted with GIRAEEE 
(WGIO- Recio-Blanco et al., in prep.), and another one devoted 
to warmer stars observed with both UVES and GIRAEEE setups 
(WG13- Blomme et al., in prep.). 

The atmospheric parameters obtained by the different WGs 
were then combined by a WG dedicated to the homogenisation 
(WG 15) that provides a final set of recommended parameters 
(Hourihane et al., in prep). More specifically, the homogenisa¬ 
tion process compares the results of the different WG for com¬ 
mon targets and calibrators and, based on that, derives and ap¬ 
plies offsets to obtain a common scale for all the Gaia-ESO Sur¬ 
vey results. We mention that the idr2/3 WG 15 homogenisation 
process has not affected the recommended values of effective 
temperature (Tgff), surface gravity (log g) and microturbolent ve¬ 
locity (^), while it has derived and applied offsets in metallicity 
[Ee/H] and in radial velocity (RV). Also, so far, the homogeni¬ 
sation flow of idr2/3 does not involve the homogenisation of the 
elemental abundances through different WG. Therefore, the ele¬ 
mental abundances used in the present work are those provided 
as recommended values by WGl 1 (Smiljanic et al. 2014). 

All the targets observed in Be 81 field have recommended 
atmospheric parameters and radial velocities in the final table of 
WG 15. The UVES targets also have element abundance de¬ 
terminations from WGll. In Table 1 we report the information 
about the UVES stars: namely, we provide their identification 
names (Cname)^ and m^, equatorial coordinates, radial veloci¬ 
ties, stellar parameters with associated errors, RV from Hayes 
& Eriel (2014) with associated errors, and membership flag (de¬ 
scribed in the Section 2.2). The idr2/3 RVs are in good agree¬ 
ment with those of Hayes & Eriel (2014). In Table 2 we show 
the information about confirmed members (based on their RV) 
observed with GIRAEEE and analysed by WGIO (see Section 
2.2 for the membership probability definition). Their stellar pa¬ 
rameters are obtained with spectral indices (Damiani et al. 2014) 
and thus the microturbulent velocity ^ is not determined. Stellar 
parameters of WG 13 are not used here because of their large 
uncertainties, and therefore they are not presented in Table 2. 

2.1. The field comparison sample 

The Solar neighbourhood Eor the MW field stars, used as 
comparison in the present work, we considered the analysis of 
WG 11. Specifically, the UVES Solar neighbourhood sample 
aims at obtaining an unbiased sample of ~5000 G-stars within 
~2 kpc from the Sun in the turn-off phase (see Gilmore et al. in 
prep, for a description of the samples). The final purpose is to 
quantify the local elemental abundance distribution functions in 
detail. Eor this paper we used the whole sample of 861 Milky 
Way solar neighbourhood stars with recommended parameters. 

The bulge and Inner disc. UVES stars observed in parallel 
with the GIRAEEE targets dedicated to the study of the Galactic 
Bulge are expected to be evolved stars belonging to both bulge 

^ the Cname is an ID used in the Gaia-ESO Survey to univocally iden¬ 
tify an object with its coordinates 

^ ID of the original photometric catalogue of Donati et al. (2014a), 
available also the CDS 
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Fig. 1. Upper left panel: in grey the Be 81 field of view of the pho¬ 
tometric data in Donati et al. (2014a), obtained with LBC@LBT; in 
black the stars observed by the Gaia-ESO survey, all inside a radius of 
6 arcmin. Lower left panel: distribution of the radial velocities versus 
distance from the cluster centre. The red line is the cluster systemic 
radial velocity. Right panel: CMD of the inner 6 arcmin, with stars ob¬ 
served indicated by larger symbols (orange squares HR09B; red crosses 
HR15N; blue circles UVES). 


and inner-disc populations. During the first 18 months of the 
Gaia-ESO Survey, 80 stars of the bulge/inner-disc were observed 
and are included in the present discussion. 

2.2. Membership 

The determination of the RVs for GIRAFFE and UVES spectra, 
respectively, can be found in Koposov et al. (in preparation) 
and (Sacco et al. 2014). We have used the recommended radial 
velocity of the UVES and GIRAFFE sample to define the cluster 
systemic velocity. 

The resulting RV distribution, including both GIRAFFE and 
UVES targets, is shown in Fig. 2. A gaussian fit of the peak 
of the distribution, identifying the cluster signature with respect 
to the field stars, gives a central velocity of 45.8 km s“' and a 
Full Width Half Maximum (FWHM) of 3.1 km s"'. To avoid a 
large contamination by field stars, we define a star as RV clus¬ 
ter member if it has a RV within l-cr from the average cluster 
radial velocity. By doing so we identify 56 RV members, 14 ob¬ 
served with UVES and 42 targeted with GIRAFFE (5 observed 
with both HR09B and HR15N, 11 with HR09B, 9 with HR15N 
and analysed by WG13, and 18 others observed with HR15N 
and analysed by WGIO). In summary, 20% of the observed tar¬ 
gets (almost 100% for UVES and 16% for GIRAFFE) have a 
RV consistent with membership. For UVES, the contamination 
by held stars has been alleviated by observing high probability 
member stars whose radial velocity were obtained by Hayes & 
Friel (2014) to assist Gaia-ESO observations. For GIRAFFE, 
the low fraction of recovered cluster members is due to the lo¬ 
cation of the cluster towards the Galactic Center where a large 
contamination by held stars is expected. 


In Fig. 3 we show [Fe/H] versus atmospheric parameters 
of the UVES stars. We remind the reader that stars observed 
with UVES are selected to be red giant stars in the clump (RC). 
They were, thus, expected to have similar atmospheric param¬ 
eters. The average error (as given in the recommended value 
table) is plotted on the left side of each panel. The l-cr and 2-cr 
around the median value areas are also marked. We note that 
the star with the lowest metallicity (19015978-0028183) is be¬ 
yond l-cr from the median Tefr and 2-cr from the median log g, 
likely being a non-member with radial velocity consistent with 
the cluster systemic velocity. This is not surprising given the 
expected high contamination of field stars. The most metal rich 
star (19013651-0027021) is also somewhat deviant and might be 
a non-member star as well. In the next sections we discuss the 
chemical homogeneity of cluster member stars and how it can 
help to disentangle member versus non-member stars. 

In Fig. 4 we plot the histogram of the Gaia-ESO recom¬ 
mended metallicities, [Fe/H], of the RV member stars observed 
with UVES and GIRAFFE. The two distributions are compat¬ 
ible: the median of the UVES stars is [Fe/H]=0.23+0.08 dex, 
whereas for the GIRAFFE stars is 0.16+0.44 dex. The UVES 
metallicities are much more concentrated in the bin centred at 
0.25 dex, while the GIRAFFE metallicities are more dispersed 
due to the lower resolution and smaller spectral range. Note that 
the errors in Table 2 are based on a single determination of stel¬ 
lar parameters, while the errors in Table 1 are computed with the 
MAD of typically ten determinations. Thus the errors on stel¬ 
lar parameters and metallicities of the GIRAFFE stars are likely 
underestimated. 

In Fig. 5 we present the HR diagram of the candidate mem¬ 
ber stars based on RV with a PARSEC isochrone (Bressan et al. 
2012) suitable for the age and metallicity of Be 81 (see next Sec¬ 
tion). We show with different colours the UVES and GIRAFFE 
stars. The most metal poor UVES star is located in the upper 
part of the plot, quite far away from the location of the other 
clump stars. The most metal rich star is located just below the 
main locus of the clump stars. The GIRAFFE targets were se¬ 
lected among main sequence (MS) stars. Several members fol¬ 
low, within the errors, the PARSEC isochrone, while other stars 
have very uncertain parameters and are not in agreement with 
the isochrone. 


3. The fundamental parameters of Be 81 

3.1. The chemical properties of Be 81 

The high-resolution and high-quality UVES spectra allow WG 
11 to determine the detailed composition of the stars of Be 81. In 
Table 3 we report the elemental abundances of the UVES stars. 
These are graphically presented in Fig. 6, where we show [El/Fe] 
of each individual UVES star. Their median values, computed 
excluding the metal-poor non-member star are shown in Table 4. 
The abundances are scaled consistently to the Solar abundances 
obtained in the analysis of WG 11. The adopted iron abun¬ 
dances are those listed in the [Fe/H] column of idr2/3 table of 
recommended parameters^. In Table 5 we show the atmospheric 
parameters and elemental abundances obtained for the selected 


cr is the median absolute deviation (MAD), computed with the Ro¬ 
bust sigma routine of IDL 

^ [Fe I/H] and [Fe II/H] abundances are available as well in the idr2/3 
table, but not for all stars of each cluster. 
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Table 1. Parameters and their uncertanities for the UVES targets in Be 81. 


Cname 

ID 

R.A. Dec. 

J2000.0 

RV 

(km/s) 

^RV 

(km/s) 

Teff 

(K) 

■^TefT 

(K) 

logg 

<^log g 

[Fe/H] 

^[Fe/H] 



RVhf14 

(km/s) 

^RVhfm 

(km/s) 

Mem. 

19013537-0028186 

32696 

285.397384 

-0.471833 

46.66 

0.40 

4966 

88 

3.09 

0.17 

0.3 

0.1 

1.49 

0.12 

- 

- 

m 

19013631-0027447 

47682 

285.401294 

-0.462442 

46.94 

0.40 

4911 

105 

2.98 

0.13 

0.23 

0.13 

1.45 

0.13 

47.1 

0.7 

m 

19013651-0027021 

33913 

285.402107 

-0.450607 

48.0 

0.40 

5082 

93 

3.42 

0.36 

0.43 

0.14 

1.4 

0.16 

49.1 

1.6 

m? 

19013910-0027114 

33770 

285.412921 

-0.453194 

48.12 

0.40 

5015 

42 

3.04 

0.12 

0.15 

0.04 

1.48 

0.09 

48.2 

0.8 

m 

19013997-0028213 

47704 

285.416526 

-0.472581 

47.73 

0.40 

4880 

73 

2.66 

0.2 

0.23 

0.22 

1.48 

0.07 

47.5 

1.2 

m 

19014004-0028129 

47688 

285.416856 

-0.470245 

48.0 

0.40 

4963 

78 

3.08 

0.41 

0.25 

0.1 

1.45 

0.19 

- 

- 

m 

19014127-0026444 

46427 

285.421947 

-0.44569 

46.28 

0.40 

4938 

47 

3.05 

0.07 

0.33 

0.11 

1.53 

0.17 

46.8 

1.1 

m 

19014194-0028172 

47703 

285.424747 

-0.471459 

47.52 

0.40 

4929 

74 

2.86 

0.14 

0.21 

0.07 

1.54 

0.11 

- 

- 

m 

19014228-0027388 

46358 

285.426169 

-0.460778 

47.72 

0.40 

4907 

57 

2.8 

0.15 

0.22 

0.08 

1.53 

0.12 

47.2 

0.7 

m 

19014498-0027496 

46477 

285.437428 

-0.463783 

47.41 

0.40 

4769 

70 

2.82 

0.18 

0.22 

0.08 

1.47 

0.12 

47.8 

0.7 

m 

19014525-0023580 

46431 

285.438563 

-0.399437 

48.73 

0.40 

5022 

12 

3.04 

0.09 

0.2 

0.08 

1.48 

0.09 

48.1 

1.3 

m 

19014769-0025108 

47674 

285.448698 

-0.419671 

47.37 

0.40 

4827 

58 

2.81 

0.26 

0.14 

0.06 

1.48 

0.19 

47.2 

0.6 

m 

19015261-0025318 

35423 

285.469182 

-0.425509 

47.85 

0.40 

5015 

83 

3.14 

0.18 

0.29 

0.14 

1.46 

0.11 

48.2 

1.0 

m 

19015978-0028183 

75596 

285.499139 

-0.471749 

46.38 

0.40 

4771 

38 

2.44 

0.08 

-0.33 

0.07 

1.45 

0.09 

46.6 

1.3 

nm 


Table 2. Parameters and their uncertainties for the GIRAFFE RV members Be 81 from WG 10. 


Cname 

Object 

R.A. Dec. 

J2000.0 

RV 

(km/s) 

^RV 

(km/s) 

SETUP 

Teff 

(K) 

^Teir 

(K) 

logg 

<5iog g 

[Fe/H] 

<5[Fe/HJ 

19012975-0027170 

42246 

285.373896 

-0.454711 

48.79 

1.39 

HR15N 

6811 

174 

3.24 

0.33 

-0.38 

0.18 

19013178-0028535 

32235 

285.382405 

-0.481533 

45.63 

0.41 

HR15N 

7253 

150 

4.09 

- 

- 

- 

19013240-0023593 

43014 

285.385012 

-0.399756 

47.83 

0.87 

HR15N 

6458 

150 

4.28 

- 

- 

- 

19013287-0031312 

41319 

285.386942 

-0.525359 

44.74 

0.51 

HR15N 

6751 

150 

4.21 

- 

- 

- 

19013372-0025163 

24010 

285.390467 

-0.421222 

46.22 

3.12 

HR15N 

6628 

142 

3.97 

0.28 

0.00 

0.12 

19013495-0028320 

32504 

285.395616 

-0.475558 

45.51 

0.72 

HR15N 

6637 

150 

4.24 

- 

- 

- 

19013809-0027105 

33782 

285.408736 

-0.452932 

46.92 

0.79 

HR15N 

6751 

150 

4.21 

- 

- 

- 

19013886-0026309 

34394 

285.411909 

-0.441952 

48.29 

3.29 

HR15N 

6545 

184 

4.35 

0.38 

0.10 

0.14 

19014212-0025400 

35245 

285.425515 

-0.427794 

48.0 

2.24 

HR15N 

6840 

126 

4.36 

0.24 

0.32 

0.11 

19014417-0025115 

45739 

285.434035 

-0.419915 

47.76 

2.2 

HR15N 

6618 

142 

4.87 

0.31 

0.23 

0.11 

19014443-0026202 

42487 

285.43514 

-0.438959 

46.77 

0.81 

HR15N 

6562 

91 

4.23 

0.21 

0.66 

0.07 

19014532-0026401 

44710 

285.438867 

-0.444473 

47.02 

0.49 

HR15N 

6637 

150 

4.24 

- 

- 

- 

19014637-0024287 

24366 

285.443198 

-0.407997 

45.48 

0.61 

HR15N 

6751 

150 

4.21 

- 

- 

- 

19014728-0025559 

35009 

285.447001 

-0.432221 

42.74 

2.32 

HR15N 

6605 

167 

5.44 

0.40 

-0.17 

0.15 

19014772-0027550 

22999 

285.448838 

-0.465302 

46.76 

1.14 

HR15N 

6011 

212 

4.65 

0.56 

-0.44 

0.18 

19014982-0028078 

22921 

285.45759 

-0.468819 

43.66 

3.45 

HR15N 

6528 

160 

4.7 

0.35 

0.16 

0.12 

19015195-0023078 

45059 

285.466451 

-0.385517 

44.37 

0.92 

HR15N 

6676 

90 

4.63 

0.2 

0.60 

0.08 

19015300-0022239 

25196 

285.470846 

-0.373304 

44.58 

0.88 

HR15N 

6947 

150 

4.17 

- 

- 

- 



4500 4600 4700 4800 4900 5000 5100 5200 
[Teff(K)] 



2.0 2.5 3.0 3.5 4.0 

[Log g] 
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- 0.2 


1.1 1.2 1.3 1.4 1.5 1,6 

[xi(km/s)] 



Fig. 3. [Fe/H] vs. atmospheric parameters of UVES stars in Be 81. The dotted lines mark the 1-cr region, while the dashed lines indicate the 2 -<t 
area. The typical error bar on each star are shown on the left side of each panel. 


Solar spectrum. For elements without a WG 11 determination, 
we scaled with the Solar abundance of Grevesse et al. (2007). 


The abundance distribution interestingly confirms the com¬ 
pactness of the cluster in most abundance ratios, whereas the 
most metal-poor star has a different chemical pattern. This is 
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Fig. 2. Distribution of radial velocities of the stars observed in the 
field of Be 81. The solid line is a gaussian fit of the distribution of the 
member stars. 



Teff(K) 

Fig. 5. The HR diagram with the PARSEC isochrone for an age of 
0.98 Gyr and a metallicity Z=0.025 (red curve) and the recommended 
Gaia-ESO parameters of UVES stars (green squares) and of GIRAEFE 
stars (blue circles). The stars that are members on the basis of their 
radial velocity but with [Fe/H] outside l-cr from the cluster average 
metallicity are shown with empty symbols. 



|Fe/H| 


Fig. 4. Distribution of [Fe/H] of the stars observed in the field of Be 81. 
The cyan rectangles show the GIRAFFE metallicities of RV member 
stars, while the blue ones indicate the UVES metallicities (including 
star 19015978-0028183). 


Table 4. Median Elemental abundances with their dispersions and aver¬ 
age errors. The non-member star is excluded from the median. 


Element 

Abundance 

(T 

< 6 > 

[Fe/H] 

0.23 

0.08 

0.10 

[O/Fe] 

0.12 

0.11 

0.05 

[Na/Fe] 

0.30 

0.10 

0.06 

[Mg/Fe] 

0.10 

0.08 

0.03 

[Al/Fe] 

0.15 

0.10 

0.08 

[Si/Fe] 

0.09 

0.09 

0.12 

[Ca/Fe] 

-0.02 

0.05 

0.09 

[Sc/Fe] 

0.07 

0.10 

0.15 

[Ti/Fe] 

0.02 

0.08 

0.13 

[V/Fe] 

0.09 

0.01 

0.11 

[Cr/Fe] 

0.02 

0.08 

0.19 

[Ni/Fe] 

-0.03 

0.08 

0.25 

[Y/Fe] 

0.03 

0.10 

0.16 

[Eu/Fe] 

0.09 

0.10 

0.21 


demonstrated in Fig. 6, which shows the dispersions around the 
median values (dashed lines) computed with the Robust sigma 
IDL routine and indicated with cr in Table 4, as well as the typ¬ 
ical error bar (black crosses in the left part of each panel and 
indicated with < 5 >) computed as the average of the individ¬ 
ual stellar errors. The error on the parameters and abundances 
of each star are discussed in Smiljanic et al. (2014), where they 
have been estimated with the method-to-method dispersion. For 
most elements these two quantities are comparable indicating a 
high degree of chemical homogeneity of the cluster, as already 
noticed for NGC 6705, Tr 20 and NGC 4815 with idrI data 
by Magrini et al. (2014). Be 81 is slightly enhanced in some 
of the (T-elements and in one of the neutron capture elements, 
having [O/Fe], [Mg/Fe], [Si/Fe] and [Eu/Fe] >0.1 dex, while it 
has [Ti/Fe] and [Ca/Fe] close to zero. In this context, the most 
metal poor star is clearly distinguishable not only for its lower 
[Fe/H] but also for differing more than l-cr in several elemental 
abundances, such as [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe], [Si/Fe], 
[Ti/Fe], [Cr/Fe], and [Eu/Fe]. Taken altogether these chemi¬ 
cal signatures indicate a clear non-membership of 19015978- 
0028183, while the most metal rich star (19013651-0027021) 
is consistent within l-cr in all panels with the median values of 
the cluster and thus can be considered a bona fide member star. 


3.2. Age, distance and reddening of Be 81 

Thanks to accurate metallicity measurements and membership 
based on RV determinations, we are able to derive an improved 
age, distance, and average reddening for Be 81 by using the clas¬ 
sical approach of isochrone fitting. As done for the open cluster 
Trumpler 20 (see Donati et al. 2014b), we chose four different 
sets of isochrones to take into account how the results change 
by adopting different models and quantify the effect of differ¬ 
ent input physics on the derived parameters: PARSEC (Bressan 
et al. 2012), BASTI (Pietrinferni et al. 2004), Dartmouth (Dot- 
ter et al. 2008), and Victoria-Regina (VandenBerg et al. 2006). 
We adopted the B, V, and I photometric catalogue described in 
Donati et al. (2014a). The accurate spectroscopic analysis of 
the cluster metallicity from the UVES spectra clearly indicates 
a super-solar chemical mixture for Be 81, in contrast to previ¬ 
ous findings which were based only on Ca-triplet spectra. Using 
this information, we are able to select a more appropriate set of 
isochrones and thus to obtain a more accurate determination of 
the cluster evolutionary status, reddening, and distance modulus. 

The best-fitting isochrone is chosen by visual examination 
as the one which can describe at the same time the main age- 
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Fig. 6. Abundance ratios of stars in Be 81. The dotted lines indicate the median+cr and median-tr regions in both [El/Fe] and [Fe/H]. The typical 
error bar on each star are shown in the left side of each panel. 


Table 5. idr2/3 Solar parameters and abundances. 


TefF 

logg 


(K) 


km 5 * 

5777±35 

4.43±0.13 

1.04+0.16 

Element 

Abundance (idr2/3) 

Abundance (G07) 

Fe 

7.53+0.06 

7.45+0.05 

0 

- 

8.66+0.05 

Na 

6.28+0.02 

6.17+0.04 

Mg 

7.71+0.06 

7.53+0.09 

Al 

6.45+0.04 

6.37+0.06 

Si 

7.46+0.04 

7.51+0.04 

Ca 

6.29+0.03 

6.31+0.04 

Sc 

3.22+0.04 

3.17+0.10 

Ti 

4.98+0.02 

4.90+0.06 

V 

3.89+0.05 

4.00+0.02 

Cr 

5.61+0.06 

5.64+0.10 

Ni 

6.15+0.14 

6.23+0.04 

Y 

2.08+0.11 

2.21+0.02 

Eu 

- 

0.52+0.06 


Abundances expressed in 12+ log(El/H) form. 


sensitive evolutionary phases: the luminosity and colour of 
the main sequence turn-off (MSTO), red hook (RH), and red 
clump (RC), defined in Donati et al. (2014a), in two different 
photometric colours {B - V and V - I, adopting Ry - 3.1, 


AbIAv - 1.29719, and AijAy - 0.60329 taken from Cardelli 
et al. (1989) and O’Donnell (1994)). We started by using the 
metallicity resulting from the UVES spectroscopic estimate, i.e. 
[Fe/H]= h- 0.23 dex (see Section 2.2), and we converted it to Z 
taking into account the different solar abundances of the four sets 
of isochrones. The difference in luminosity between the MSTO 
and RC is adopted as the most reliable constraint on the cluster 
age. In fact, while the colour difference between the same two 
indicators is related to age as well, it is used here as secondary 
age indicator because colour properties are more affected by the¬ 
oretical uncertainties, like colour transformations and the super- 
adiabatic convection, while luminosity constraints are more reli¬ 
able. The errors on the estimated cluster parameters are mainly 
due to the uncertainties in the definition of the age indicators be¬ 
cause of the heavy contamination by field interlopers, while the 
photometric error does not have a significant impact on the er¬ 
ror budget (see Donati et al. 2014a). We make use of the RV 
to select the most probable cluster members among the many 
stars observed in order to highlight the most probable cluster se¬ 
quences in the CMD and to have more stringent constraints on 
the best-fitting models. Moreover, we tweak the age determina¬ 
tion looking for the best-fitting isochrone in the theoretical plane 
Teff, log g by using the RC star parameters of the UVES targets 
(see Table 1). 
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For the PARSEC set (see Fig. 7, left panel), for which 
Zq = 0.015, we adopted Z=0.025. The luminosity and colour 
of the RC are well reproduced as are those of the RH and upper 
MS. In Fig. 5 we show, as an example, the isochrone fit in the 
theoretical plane which, together with the photometric data, was 
used to determine the best age, i.e. 1.0+0.1 Gyr. We note that 
in the case of the V,V - I CMD the RC and MS colours of the 
model are slightly bluer than observations, and this is evident 
also in the best-fit with the three other isochrones. As mentioned 
above, this effect may be due to uncertainties in the temperature- 
colour transformations or to the differences in chemical mixtures 
between models and observations or, finally, to small systematics 
in the B and I photometry that in principle might have a larger 
impact on colours, of the order of some hundredths of magni¬ 
tude. Considering, however, that metallicity is fixed by spec¬ 
troscopy and age by the difference in luminosity between RC 
and TO, we actually used V,B - V and V,V - I photometry to 
adjust the distance modulus and reddening, and deemed this fit 
quite satisfying and the small mismatch of the RC colour accept¬ 
able. In the case of the BASTI isochrones (Z=0.03), the best-fit 
is obtained for a younger age of 0.75 + 0.1 Gyr. The Dartmouth 
isochrones and the Victoria-Regina isochrones do not include 
the evolved phases after the RGB for the age of Be 81 and thus 
the cluster’s parameters for these two models are more loosely 
constrained. We proceeded by choosing the best-fitting solution 
as the one that best matches the MS and MSTO morphologies 
(see Fig. 7, right panel). We used models with [Fe/H]=+0.23 for 
both sets of isochrones. The results obtained with the four sets 
are given in Table 6, where we indicate age, distance modulus, 
reddening, distance from the Sun and the Galactocentric distance 
from the Galactic plane, and mass of the stars at the MSTO. The 
cluster’s centre coordinates are taken from Donati et al. (2014a): 
RA=19:01:42.8 and Dec=-0:27:07.7. 

Since from spectroscopic analysis it turned out that Be 81 
stars show an a-enhanced chemical pattern ([a/Fe]—1-0.2, see 
Section 3.1), we look also at the model predictions when an a- 
enhanced mixture is taken into account. Slightly younger ages 
are obtained. The results are reported in Table 6. 

With respect to literature values (see Donati et al. 2014a; 
Sagar & Griffiths 1998), we found an acceptable agreement with 
the determination of the age, distance and reddening of the clus¬ 
ter. 

4. A comparison with inner disc ciusters and fieid 
stars 

The large variety of available chemical elements, belonging 
to different nucleosynthesis channels, allows us to compare, 
in a fully consistent way, the detailed chemical patterns of 
open clusters and of field stars observed and analysed homo¬ 
geneously by the Gaia-ESO Survey. Besides Be 81, three other 
old/intermediate-age clusters are available in idr2/3. For all of 
them, namely Tr 20, NGC 6705, and NGC 4815, the stellar pa¬ 
rameters, abundances, and radial velocities, were also provided 
in the first data release, idrI, and analysed in several papers: 
Donati et al. (2014b) derived membership, age, distance, and 
reddening of Tr 20; Friel et al. (2014) derived the parameters of 
NGC 4815 and studied its chemical composition; Cantat-Gaudin 
et al. (2014) studied the properties of NGC 6705 and searched 
for possible inhomogeneities in its chemical composition; finally 
Magrini et al. (2014) searched for unique signatures in the chem¬ 
ical patterns of the three clusters as possible indication of their 


birthplace. Their CNO abundances have been discussed by Taut- 
vaisiene et al. (2015) using idr2/3 results. Their main parameters 
are summarised in Table 7. 

Here we use the idr2/3 stellar parameters and abundances 
of NGC 4815, NGC 6705, and Tr 20 to compare their chemi¬ 
cal patterns with Be 81 and with field stars. Although shown in 
Fig. 6, in the following we do not discuss Al and Na in terms 
of Galactic chemical evolution. They are indeed potentially in¬ 
teresting elements for comparison with the field samples and de¬ 
duction of the chemical evolution of the inner disc. However, the 
abundances derived in idr2/3 are all computed in local thermo¬ 
dynamic equilibrium (FTE), and those elements are affected by 
relevant non-local thermodynamic equilibrium (NFTE) effects. 
All stars observed with UVES in Be 81 are in the same evolu¬ 
tionary phase (RC), so that the internal dispersion in abundance 
is not affected by this problem. However, the absolute value of 
Na and Al abundances and the comparison with stars of different 
parameters could be unreliable without these corrections. Fur¬ 
thermore, both Na and Al are subject to (extra-)mixing processes 
in stars of the mass involved here (see, e.g., Charbonnel & Fa- 
garde 2010; Karakas & Fattanzio 2014). These elements will be 
discussed in a forthcoming paper (Smiljanic et al., in prep) and 
will not be used. 

As is well known, the ratio of a-elements to iron, [a/Fe], is 
commonly used to trace the star-formation timescale in a sys¬ 
tem, because it is sensitive to the ratio of SNII (massive stars) 
to SNIa (intermediate mass binary systems with mass transfer) 
that have occurred in the lifetime of the system. However, the 
so-called a-elements, that are generally associated with the ex¬ 
plosion of massive stars, do not have exactly the same place of 
production, and some of them, as Si, Ca, and Ti, have also a 
significant contribution from SNIa. In addition, there are also 
differences in the way these elements are produced in massive 
stars: O and Mg are produced during the hydrostatic He burn¬ 
ing, and their yields are not expected to be significantly affected 
by the SNII explosion conditions, while Si, Ca and Ti are mostly 
produced during the SNII explosion. This distinction is also seen 
in the observations (e.g., Fulbright et al. 2007), where Si, Ca and 
Ti usually track one another, but O and Mg often show different 
trends with [Fe/H]. The "pure" a-elements (O, Mg) are produced 
almost entirely by massive stars, and thus in shorter time-scales. 
In the bulge, the observed Mg overabundances and the general 
a-enhancement agree with the predictions of the chemical evo¬ 
lution models of Matteucci & Brocato (1990) and Grieco et al. 
(2012) and they are explained with a rapid formation time scale 
for the bulge. A summary of the place of production and of the 
contribution of SNII and SNIa is shown in Table 8. The Table is 
computed using the information given in Woosley et al. (2002) 
and in Iwamoto et al. (1999, model W7) and summarises the 
origin of the main isotope of each element. 

In Figure 9 we show the median [El/Fe] abundance ratios 
versus the atomic number in the four old/intermediate-age clus¬ 
ters located within 7 kpc from the Galactic Centre, whereas 
in Figure 10 we show the chemical pattern of the four clus¬ 
ters compared to the median abundance ratios of stars with 
-0.10<[Fe/H]< +0.33 dex in the solar neighbourhood [382 
stars] and in the inner-disc/bulge [42 stars] samples, i.e., embrac¬ 
ing the maximum metallicity range of the four clusters within 1- 
cr. We have excluded all stars with errors on [Fe/H] larger than 
0.15 dex. 

In both figures, the available elements include the a-elements 
(O, Mg, Si, Ca, Ti), the refractory elements (Sc, V), the iron- 
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Fig. 7. CMD obtained for stars inside 3' using the photometry from Donati et al. (2014a) and best-fit isochrone for different evolutionary models 
(from left to right: PARSEC, BASTI, DARTMOUTH, and Victoria-Regina). GES target non members are plotted as black crosses, while members 
are highlighted with orange (GIRAFFE) and blue (UVES) points. See Table 6 for the adopted parameters for the isochrone fitting. 



peak elements (Cr, Ni), and the neutron capture elements (Y, 
Eu). From Fig. 9, we can make some general considerations on 
the similarities in the chemical patterns of the inner clusters. 


Considering the errors that affect the abundance ratios, 
we define an element as enhanced with respect to Solar if 
[El/Fe]>0.1 dex. Also, we consider that the abundance of an el¬ 
ement is globally enhanced in the inner cluster sample if at least 
three of the four clusters have [El/Fe]>0.1 dex for that element. 


NGC 6705, NGC 4815, and Be 81 have enhanced abun¬ 
dances of the a-elements O, Mg, and Si, with NGC 6705 hav¬ 
ing the highest abundance ratios of [Mg/Fe] and [Si/Fe] while 
[O/Fe] does not exceed a 0.1 dex enhancement. The group of 
elements from Ca to Ni shows Solar abundance and they have 
similar values in the four clusters. Y and Eu are very similar 
in all clusters, close to Solar value for [Y/Fe] and slightly super¬ 
solar for [Eu/Fe]. From Fig. 10, we find that the chemical pattern 
of Be 81 is nearly identical to that of inner-disc/bulge stars: both 
populations show similar a-enhancements for the lower atomic 
number a elements (O, Mg, and Si), while Ca, Ti and iron-peak 
elements, are almost solar. For comparison, the Solar neighbour¬ 
hood stars in the same metallicity range are plotted; they all have 
[Mg/Fe], [Si/Fe] and [O/Fe] close to zero. 


5. The inner clusters and the disc evolution 

For the first time, thanks to the Gaia-ESO Survey, we have at our 
disposal a significant number of open clusters together with field 
stars, homogeneously analysed and for which elements from dif¬ 
ferent nucleosynthesis channels are available. They allow us to 
analyse the typical chemical patterns in different regions of the 
disc of our Galaxy. Clusters are unique tools because their dis¬ 
tances and ages can be derived with great accuracy, allowing 
reliable placement in space and time in our Galaxy. The fam¬ 
ily of clusters here described belongs to an important but poorly 
investigated region of our Galaxy: the inner disc. 

There are two main results; 

1. As noticed in Magrini et al. (2014), there are ’real’ differ¬ 
ences in the abundance patterns of the clusters even if they 
are located at similar distances from the Galactic centre. One 
very noticeable difference is the higher o'-enhancement of 
NGC 6705, that might suggest an inner birthplace for this 
cluster, or, as we see in next sections, possibly an episode of 
local enrichment. 

2. Despite the small but ’real’ differences, a common behaviour 
of the family of old/intermediate-age clusters analysed so far 
can be seen: they are almost solar in most of the iron-peak 
and heavy a-elements, while they are enhanced in the light a- 
elements O, Mg, and Si. We discuss in the next sections how 
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Fig. 8. The same as Figure 7 but in the V,V - I CMD 


much this enhancement is real or due to Non-Local Thermo¬ 
dynamical Equilibrium (NLTE) effects. 

5.1. Is the a-enhancement seen in the inner clusters real? 

Three of the clusters of our inner disc sample have [O/Ee], 
[Si/Ee], and [Mg/Ee]>0.1 dex. [Si/Ee] is very high in NGC 6705, 
while it is around 0.1 dex for the remaining two clusters. Eirst, it 
is necessary to investigate if the O, Mg, and Si enhancements are 
real or related to other effects, such as NLTE. All stars analysed 
in the four clusters are indeed giant stars, and thus NLTE might 
affect the abundance of some elements. If the NLTE effect is 
large, the comparison between the giant stars in clusters and the 
Sun or the Milky Way field solar-neighbourhood dwarf stars can 
be misleading. We clarify that our ’definition’ of NLTE is totally 
empirical; in the following we indicate as a ’NLTE effect’ the 
sum of all un-modelled physics plus potential model biases. 

Oxygen. In Be 81, as well as in the other clusters and field 
stars, oxygen abundance was determined from the forbidden 
[O I] line at 6300.31 A, as described in detail by Tautvaisiene 
et al. (2015). Lines of [O I] are considered as very good indica¬ 
tors of oxygen abundances because they are not only insensitive 
to NLTE effects, but also give similar oxygen abundance results 
with 3D and ID model atmospheres (cf. (Asplund et al. 2004); 
(Pereira et al. 2009)). Thus, the comparison of giants and dwarfs 
should not be affected by any such bias. 

Magnesium. Mg is less affected than the other a-elements 
by model dependent differences in the yields of SN II with pro¬ 


genitors of different masses. Therefore it constitutes one of the 
best tracers of the evolutionary timescales of the Galaxy. How¬ 
ever, the NLTE effects in Mg are not easily quantifiable. From 
the theoretical point of view. Merle et al. (2011) computed the 
Mg I abundance corrections for several lines in the optical in¬ 
terval as functions of atmospheric parameters for cool evolved 
stars. Their NLTE corrections are computed as a function of 
stellar parameters and metallicity. At solar metallicity for a star 
with Tefr=5000 K and log g=2 the NLTE corrections of differ¬ 
ent Mg I lines vary from -0.03 to 0.18 dex. From the observa¬ 
tional point of view, several works in the literature have anal¬ 
ysed dwarf and giant stars in the same cluster to investigate if 
there is any difference in the abundance ratios related to the type 
of the star analysed. A study of dwarf and giant stars in the 
Hyades open cluster (Schuler et al. 2009) showed that Mg abun¬ 
dances are overabundant by 0.2-0.5 dex in the giants relative to 
the dwarfs. Such large enhancements are not predicted by stellar 
models, and they explained them as likely due to NLTE effects, 
significantly larger than those available in the literature. A sim¬ 
ilar result was obtained by Yong et al. (2005) comparing abun¬ 
dances in dwarf and giant stars in M 67. However, Pasquini et 
al. (2004) and by Blanco-Cuaresma et al. (2015) found contrary 
results for the open cluster IC 4651 ([Fe/H] = 0.11). 

Silicon. The few studies on its NLTE behaviour indicate 
strong effects only in low metallicity stars. Recently, Bergemann 
et al. (2013) found that the effect of NLTE is a reduction of the 
number densities of atoms in a minority ionisation stage, thus, 
for a given abundance, the line equivalent width is smaller than 
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Fig. 9. The abundance ratios of the idr2/3 old and intermediate-age open clusters located in the inner disc as a function of the atomic number. In 
black are the abundance ratios of Be 81, in yellow NGC 6705, in green NGC 4815, and in red Tr 20. The dashed line marks the solar abundance, 
and the shaded area the l-cr regions, with cr=0.1 dex. 


that given by LTE. This effect is larger in metal-poor stars. At 
solar metallicity, the opposite behaviour is encountered, show¬ 
ing negative NLTE effects such that the LTE abundance is higher 
than the NLTE one. 

In order to better investigate the issue for Si and Mg, we test 
the nature of the a-enhancement in our results using the results 
for one of the Gaia-ESO calibrator open clusters (see Pancino 
et al. in prep.), M 67. The data for M 67 are, indeed, perfectly 
suited because they contain stars belonging to different evolu¬ 
tionary phases that have been analysed with the same tools and 
methods as the inner-disc clusters and field stars of idr2/3. 

In Eig. 11, we show [Mg/Ee] and [Si/Ee] vs Teu in mem¬ 
ber stars of M 67 analysed in idr2/3. The +l-cr and +2-cr re¬ 
gions from the average are shown. The sample of stars with 
Teff >5700 K includes TO, MS, and SG stars. The stars with 
Teff <5700 K are instead giant stars with log g ranging from 
2.97 to 3.46. We note that for both elements a slight overabun¬ 
dance in giants is present with respect to TO, MS, and SG stars. 
The average [Mg/H] for stars with Teu <5700 K is 0.04+0.02 
and for [Si/H] is 0.08+0.01, while the stars with Teff >5700 K 
have average values of -0.09+0.04 and -0.03+0.03 from [Mg/H] 
and [Si/H], respectively. We can use these differences to roughly 
estimate the amount of NLTE effects: 0.13 dex and 0.11 dex for 
[Mg/H] and [Si/H], respectively. If we consider the abundances 
with respect to iron, the amount of NLTE effects is 0.12 dex and 
0.09 dex for [Mg/Ee] and [Si/Ee], respectively. 

It is worth mentioning that M 67 has been used by One- 
hag et al. (2014) to study the effect of selective atomic diffusion 
of chemical elements in stars in different evolutionary phases. 
More specifically, they investigated the detailed chemical com¬ 


position of 14 stars located on the main sequence, at the turn-off 
point, and on the early sub-giant branch, finding that the heavy- 
element abundances are reduced in the TO and main sequence, 
by typically <0.05 dex, when compared to the abundances of the 
sub-giants (even lower for ratios of metals). We note that the 
effect of diffusion in our data is circumscribed in the +lcr area. 

We finally mention that Blanco-Cuaresma et al. (2015) study 
the differences in the abundance ratios in cluster giant and 
dwarf member stars. They show the chemical differences be¬ 
tween dwarfs and giants in several clusters, including IC 4651, 
M 67, NGC 2447, and NGC 3680. They find that, on average, 
[Mg/Ee] and [Si/Ee] are enhanced in giant stars by ~0.06 dex and 
~0.10 dex with respect to dwarfs. In particular, for M 67, the ef¬ 
fect they measured is 0.07 and 0.10 dex in [Mg/Ee] and [Si/Ee], 
respectively, thus comparable with what we find with Gaia-ESO 
recommended values. 

In conclusion, for oxygen, no NLTE has to be applied, and 
the three clusters for which this element has been measured show 
a [O/Ee]>0.1 dex. The enhancement in [Mg/Ee] and [Si/Ee] ob¬ 
served in the inner-disc/bulge population and in Be 81, Tr 20, 
and NGC 4815 with respect to the Solar neighbourhood dwarf 
population can be related to different evolutionary phases of the 
observed stars (see Eigs.9 and 10). On the other hand, [Mg/Ee] 
and [Si/Ee] are remarkably higher in NGC 6705 than in the other 
inner-disc clusters and in the field populations. We can evidently 
exclude that this enhancement is due solely to NLTE effects. The 
abundance ratios of NGC 6705 clearly indicate an a-enhanced 
population whose origin needs to be investigated. 
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Fig. 10. The abundance ratios vs atomic numbers of the idr2/3 stars in the solar neighbourhood and in the inner disc/bulge (selected in the 
[Fe/H] range from -0.10 to -fO. 33 dex, i.e., including the metallicity ranges of the four clusters) compared with the chemical pattern of NGC6705, 
Trumpler 20, NGC 4815, and Be 81. The continuous lines mark the Solar abundances, and the shaded area the l-cr region, with cr=0.1 dex. 
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Fig. 11. The abundance ratios [Si/H] and [Mg/H] vs T^ff of member stars in M 67 from idr2/3 analysis. The continuous lines mark the ±i-a 
from the average value, whereas the dotted lines indicate the ±2-cr region. 


5.2. Comparing the inner dusters with high-a metai-rich 
popuiation 

Once shown that the enhancements of [O/Fe], [Si/Fe] and 
[Mg/Fe] in NGC 6705 are genuine, we can try to compare its 
chemical pattern with those of a population that presents simi¬ 
lar chemical characteristics, the high-a metal-rich stars (HAMR) 
with [Fe/H]=0.03+0.12. 


The HAMR stars-whose a-enhancement, computed with 
Mg, Si, and Ti, is in the range 0-0.2 dex- are rare objects in 
the solar neighbourhood, and their exact nature is still under de¬ 
bate. Adibekyan et al. (2011) identified these stars and suggested 
that they may form a distinct family, based on a “gap” in metal¬ 
licity between HAMR and classical thick disc stars. Other au¬ 
thors also confirmed the existence of the a-enhanced metal-rich 
stars but with different conclusions on their origin (e.g. Gazzano 
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et al. 2013; Nissen et al. 2014; Bensby et al. 2014; Adibekyan 
et al. 2014). A comparison of the main physical properties of 
HAMR stars with those from the simulations of Roskar et al. 
(2012) suggests that the HAMR stars might have their origin in 
the central part of the Galaxy, within about 2 kpc (Adibekyan 
et al. 2013). Looking at Figs. 9 and 10, one would claim that 
stars in Be 81 and in the inner-disc clusters are similar to the 
a-enhanced metal-rich stars discovered in the solar neighbour¬ 
hood. They have indeed solar or super solar [Fe/H] (see Ma- 
grini et al. 2014, and Table 4 for Be 81) and they are, in general, 
enhanced in the light a-elements, as O, Mg, and Si. At a first 
glance, the chemical patterns of Be 81 shown in Fig. 9 pretty 
much resemble the metal rich a-enhanced population identified 
in the solar neighbourhood. Note that HAMR stars studied by 
Adibekyan et al. (2011) have Tgff > 4900 K and only 5 stars 
have log g<3.8 dex, therefore their sample is essentially com¬ 
posed of dwarf stars. If confirmed, the existence of HAMR stars 
in an open cluster at 5-7 kpc from the Galactic Centre, might put 
important constraints on the stellar migration and in the birth¬ 
place of HAMR stars. 

However, we should consider that there are important differ¬ 
ences between HAMR stars and the inner disc clusters. First, 
there is a large difference in age between the two populations. 
The age of HAMR stars in the solar neighbourhood is not pre¬ 
cisely determined. However, they are on average believed to be 
quite old. Adibekyan et al. (2011) showed that the HAMR and 
thick disc family stars have almost the same age, being on av¬ 
erage older than thin disc stars by about 3 Gyr. At the same 
time, they found that HAMR stars are quite different from the 
thick disc stars being more metal-rich and with different chemi¬ 
cal patterns. They also show properties in common with metal- 
rich bulge stars (e.g., Fulbright et al. 2007; Bensby et al. 2014; 
Ness et al. 2013). On the other hand, the age and distance of 
NGC 6705 are well known by isochrone fitting (see Table 7), 
and it is much younger. 

Second, in a ’classical’ chemical evolution view of the Milky 
Way disc (e.g. Tosi 1988; Matteucci & Francois 1989; Chiap- 
pini et al. 1997, 2000; Portinari et al. 1998; Boissier & Prantzos 
1999) we do not expect that relatively young stars born in the in¬ 
ner disc to show an a-enhancement. It is instructive to use a re¬ 
cent and simple chemical evolution model that does not include 
the effect of radial migration to study the most likely birthplace 
of the HAMR stars. In Fig. 12 we show a plot of [Fe/H] and 
[a/Fe] vs. Rgc in the chemical evolution model of Magrini et 
al. (2009). The three curves show the prediction for [Fe/H] and 
[a/Fe] as a function of Rgc at present, 5 Gyr ago and 10 Gyr ago. 
The model predicts that the formation of a family of a-enhanced 
metal-rich stars is very limited in space and in time. They formed 
in the inner part of the disc and during the first epochs of disc for¬ 
mation (black curve), then the contribution of SNIa decreased 
the [a/Fe] ratio to solar values. Thus, further stellar populations 
have higher [Fe/H] but lower [a/Fe] than the oldest ones. Similar 
results can be obtained even with more sophisticated models, as 
those by Minchev et al. (2013), where the combination of high 
[a/Fe] and young ages is not predicted. For instance, in Figure 2 
of Minchev et al. (2013), stars with [O/Fe]=0.2 are older than 
7 Gyr, whatever their birth location. Stars with slightly smaller 
[O/Fe] can be younger at the condition that they are born in the 
very outer disc, but then they would be metal poor. 

We plot also [Fe/H] and [a/Fe] for the four clusters and 
for the solar neighbourhood and inner-disc samples in the same 
metallicity range of the clusters. We compute [a/Fe] consider¬ 


ing [O/Fe], [Mg/Fe] and [Si/Fe]-the latter two corrected for the 
NLTE effects. For the field populations we assumed a distance 
of 8+1 kpc for the Solar neighbourhood stars, and of 5+2 kpc 
for the inner-disc sample. The uncertainties on the x-axes reflect 
their interval of distances. 

The four clusters are in good agreement with the model pre¬ 
diction in the [Fe/H] vs. Rgc plane, however in the [a/Fe] vs. 
Rgc plane NGC 6705 is displaced with respect to what is ex¬ 
pected for its age. Stars with the age of ~0.3 Gyr, indeed, should 
have solar or sub-solar [a/Fe] given their age and place of birth 
(blue curve of Fig. 12) and thus their a-enhancement is totally 
unexpected. Also the location of Be 81 is a bit surprising: as 
we discussed above, this cluster does not have a high [a/Fe] ‘per 
se’. However, if we consider its [a/Fe] in conjunction with its 
[Fe/H], we find it a bit higher than expected. On the other hand, 
[a/Fe] in Tr 20 and NGC 4815 are consistent, within the errors, 
with the model curves (in between red and blue curve of Fig. 12). 

Therefore, the stars in NGC 4815 and Tr 20 have a [a/Fe] 
that is fully consistent with their age and location in the disc as 
indicated by the predictions of the chemical evolution model of 
Magrini et al. (2009). Be 81 has a slightly higher [a/Fe] than ex¬ 
pected for its [Fe/H] and Rgc, However, the most striking case is 
that of NGC 6705, which shows a genuine enhancement [a/Fe] 
unlikely related to the HAMR stellar population in the solar 
neighbourhood, since they were born in a completely different 
epoch. Due to the young age of NGC 6705, its a-enhancement 
cannot be simply explained by disc chemical evolution, and also 
the hypothesis of Magrini et al. (2014) of radial migration from 
the inner disc is unlikely due to the young age of this cluster. 

Very recently, several studies revealed the presence of rel¬ 
atively young a-enhanced stars. Chiappini et al. (2015) re¬ 
port the discovery of a group of apparently young CoRoT red- 
giant stars exhibiting enhanced [a/Fe] abundance ratios. These 
stars are more numerous in their inner-disc sample. Also Mar- 
tig et al. (2014) analysing a sample of giants with seismic pa¬ 
rameters from Kepler found a small percentage of stars with 
[a/Fe] >0.2 dex and ages below 4 Gyr. As in the case of our 
inner disc clusters, the existence of this kind of star cannot be 
explained in the framework of standard chemical evolution mod¬ 
els of the Milky Way. Their presence might indicate that the 
chemical-enrichment history of the Galactic disc is more com¬ 
plex. Chiappini et al. (2015) analyse some possible interpreta¬ 
tions of the nature of these stars: i) they might have been formed 
close to the end of the Galactic bar, near co-rotation where gas 
stays inert for longer times and where the mass return from older 
inner-disc stellar generations is expected to be highest with the 
net effect of in-situ gas dilution; ii) they might originate from a 
recent gas-accretion event. However, there is a substantial differ¬ 
ence between the stars discussed in Chiappini et al. (2015) and 
the inner-disc open clusters: Chiappini’s stars tend to be more 
metal poor than our clusters, with [Fe/H] ranging from ~ -0.6 
to ~0.1 dex, and they are on average older (~4 Gyr). Might the 
inner-disc open cluster be a younger tail of the "young’ metal 
poor a-enhanced population? Is the high [a/Fe] of both popula¬ 
tions related to their place of birth, near to co-rotation? 

In the next section, we try to estimate another possibility: 
the effect of a very local enrichment in the molecular cloud from 
which NGC 6705 was bom to explain its chemical pattern. 
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Table 9. Cluster enrichment by SN II using yields of WW95. 


Mass 

(Mo) 

[O/Fe] 

[Mg/Fe] 

[Si/Fe] 

15 

0.09 

0.03 

0.23 

18 

0.16 

0.09 

0.27 

25 

0.39 

0.19 

0.49 


5.3. The local enrichment hypothesis: the effect of a 
supernova explosion 

Can a single SNII explosion explain the chemical pattern of 
NGC 6705 with its high [a/Fe]? Using Gaia-ESO data, Berge- 
mann et al. (2014) found a large scatter of [Fe/H] and [Mg/Fe] 
of stars at each age towards the inner disc. Malinie et al. (1993) 
explained the observed scatter at each Galactocentric radius and 
age as due to the fact that the disc contains chemical inhomo¬ 
geneities that can avoid remixing long enough. This allows the 
activation of the star formation process in both enriched and un¬ 
enriched portions, with typical timescales of 10*^-10^ yr. Gilmore 
et al. (1990) suggested that the scatter could be a consequence of 
self-enrichment in giant molecular cloud complexes. Their life¬ 
time is short compared to the time scale of Galactic evolution, 
but much larger than the lifetime of the massive stars. Therefore, 
self-enrichment could occur inside a giant molecular cloud, and 
after the disruption of the cloud the metallicity inhomogeneity 
survives. It spreads only through a Galactic annulus, but it can 
be easily detected if the cluster formed from the molecular cloud 
is still observable. 

Starting from the assumption that the interstellar medium 
from which NGC 6705 was bom had a mass comparable to its 
present time mass (~3500 Mq, Cantat-Gaudin et al. 2014), we 
compute the enrichment due to the explosion of a massive star 
of 15, 18 or 25 Mq (yields from Woosley & Weaver (1995)) in 
the hypothesis of a uniform and instantaneous recycling of the 
medium. The enrichment has been computed assuming Solar 
initial composition (Grevesse et al. 2007). 

In Table 9, we present the results of this calculation; the ex¬ 
plosion of a massive star with a mass 15-18 Mq might be able 
to explain the local enrichment of NGC 6705. The yields of 
Woosley & Weaver (1995) can also explain the higher enrich¬ 
ment in [Si/Fe] with respect to the other elements. The hy¬ 
pothesis of local enrichment is supported also by Roy & Kunth 
(1995) where characteristic timescales for homogenisation are 
short, 10^-10*^ yr on scales of 100-1000 pc. The number of stars 
with a mass in the range 15-18 Mq computed with the initial 
mass function of Cantat-Gaudin et al. (2014) in the whole area of 
the cluster is about 2. Given the present total mass of NGC 6705, 
the expected number of stars with 15 Mq <M< 18 Mq is thus 
consistent with an early enrichment by one or more type II super¬ 
novae with such massive progenitors. However, the hypothesis 
that the SN was a cluster member implies at least three assump¬ 
tions; namely, i) a prolonged star formation (at least 10^ Myr : 
timescales of SNII with progenitors of 15 Mq <M< 18 Mq), 
where massive stars have formed first and have enriched the 
medium from which the lower mass stars have subsequently 
formed; ii) that the explosion of the SN has not dissipated the gas 
from which the low mass stars have formed; Hi) that no low mass 
stars have formed from the pristine (not yet enriched) gas, or that 
those stars have escaped the clusters. Similar problems are faced 
also in globular cluster formation and in their self-enrichment. 


6. Summary and conclusions 

We analyse the Gaia-ESO results (idr2/3) for four inner disc 
clusters, including the old cluster Be 81 which has not been stud¬ 
ied in the Gaia-ESO survey so far. Using the improvements ob¬ 
tained in the metallicity measurement and the membership with 
RV determinations, we derive age, distance, and average red¬ 
dening of Be 81 by using the classical approach of isochrone fit¬ 
ting. The best-fit isochrone gives an age of ~ 1 Gyr, a reddening 
E(B - V) =0.83 and a Galactocentric distance of 5.45 kpc. We 
compare the chemical pattern of Be 81 and of the other three old 
open clusters of idr2/3 with the solar-neighbourhood and inner- 
disc/bulge stars. We find that the inner-disc/bulge giant stars and 
three open clusters show an enhancement in the light a-elements 
(O, Si, Mg). We use the calibrator cluster M 67 to infer how 
much of this enhancement is related to NFTE effects, finding 
that the enhancement of [Mg/Fe] and [Si/Fe] in giants with re¬ 
spect to dwarfs is ~0.10 dex. Applying these empirical NFTE 
corrections to the chemical patterns of open clusters, we find 
that Tr 20, NGC 4815, and Be 81 do not show any remarkable 
a-enhancement. On the other hand, NGC 6705 still shows an un¬ 
common pattern, with a significant a-enhancement. We compare 
the cluster [a/Fe] with those of the HAMR stars, finding that they 
belong to very different epochs of the lifetime of the Galaxy. We 
relate them also to the recently discovered population of young 
(age~4 Gyr) and a-enhanced stars ((Martig et al. 2014), (Chi- 
appini et al. 2015)) mainly located towards the Galactic Centre. 
The inner-disc clusters, with ages between ~0.3 to ~1.5 Gyr, 
and the "young" a-enhanced stars might belong to an evolution¬ 
ary sequence of the same population, born in the inner part of the 
disc. The common chemical features of some of the clusters with 
the "young" a-enhanced stars might indicate peculiar conditions 
for the star formation near to the co-rotation. In our discussion, 
we also tentatively explain the very high [O/Fe], [Mg/Fe], and 
[Si/Fe] in NGC 6705 with an episode of local enrichment due 
to a type II supernova in the mass range 15-18 Mq, as, e.g., 
contemplated by Malinie et al. (1993) to explain the observed 
inhomogeneities at each Galactocentric radius. 
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Table 3. Abundances in the 12+log(X/H) form of the UVES stars in Be 81 from idr2/3 recommended values. 


Cname 

Oi 

Nat 

Mg I 

Alt 

Sit 

Cat 

Sc I 

19013537-0028186 

9.11+0.05 

6.86+0.31 

8.05+0.29 

6.82+0.09 

7.80+0.12 

6.54+0.13 

3.54+0.04 

19013631-0027447 

9.01+0.05 

6.81+0.07 

8.05+0.08 

6.83+0.02 

7.81+0.15 

6.45+0.08 

3.69+0.15 

19013651-0027021 

9.26+0.05 

6.91+0.05 

8.09+0.15 

6.92+0.03 

7.81+0.18 

6.72+0.10 

3.64+0.17 

19013910-0027114 

8.86+0.05 

6.75+0.06 

7.96+0.02 

6.77+0.06 

7.70+0.09 

6.45+0.10 

3.44+0.17 

19013997-0028213 

9.05+0.05 

6.67+0.26 

8.00+0.17 

6.71+0.19 

7.77+0.14 

6.44+0.07 

3.42+0.15 

19014004-0028129 

9.01+0.05 

6.90+0.06 

8.06+0.01 

6.83+0.06 

7.83+0.11 

6.55+0.10 

3.61+0.08 

19014127-0026444 

9.01+0.05 

6.78+0.35 

8.01+0.10 

6.81+0.03 

7.81+0.12 

6.44+0.13 

3.54+0.20 

19014194-0028172 

8.87+0.05 

6.81+0.05 

8.07+0.13 

6.83+0.01 

7.79+0.11 

6.43+0.07 

3.49+0.12 

19014228-0027388 

8.88+0.05 

6.86+0.06 

8.20+0.01 

6.84+0.01 

7.77+0.12 

6.47+0.07 

3.52+0.19 

19014498-0027496 

9.11+0.05 

6.70+0.29 

7.95+0.09 

6.77+0.02 

7.66+0.13 

6.49+0.09 

3.47+0.21 

19014525-0023580 

8.96+0.05 

6.81+0.13 

7.96+0.14 

6.83+0.09 

7.75+0.09 

6.50+0.07 

3.46+0.08 

19014769-0025108 

8.96+0.05 

6.82+0.04 

8.03+0.07 

6.83+0.03 

7.81+0.19 

6.44+0.05 

3.67+0.11 

19015261-0025318 

8.96+0.05 

6.84+0.29 

8.11+0.01 

6.78+0.08 

7.76+0.12 

6.54+0.09 

3.58+0.15 

19015978-0028183 

8.83+0.05 

6.04+0.06 

7.66+0.01 

6.37+0.07 

7.32+0.08 

5.96+0.10 

3.04+0.11 

Cname 

Til 

Vi 

Cri 

Nit 

Yii 

Eu It 


19013537-0028186 

5.27+0.18 

- 

5.99+0.14 

6.52+0.29 

2.53+0.12 

- 


19013631-0027447 

5.18+0.11 

4.21±0.11 

5.79+0.22 

6.55+0.25 

2.51+0.16 

0.73+0.26 


19013651-0027021 

5.36+0.17 

- 

6.04+0.17 

6.62+0.32 

2.50+0.22 

- 


19013910-0027114 

5.09+0.11 

4.13+0.11 

5.85+0.23 

6.34+0.21 

2.40+0.32 

0.76+0.13 


19013997-0028213 

4.99+0.16 

- 

5.71+0.23 

6.31+0.34 

2.49+0.12 

- 


19014004-0028129 

5.19+0.09 

4.25+0.11 

5.88+0.19 

6.52+0.26 

2.55+0.16 

0.93 ±0.21 


19014127-0026444 

5.16+0.13 

- 

5.87+0.16 

6.49+0.19 

2.47+0.23 

- 


19014194-0028172 

5.18+0.13 

- 

5.95+0.23 

6.50+0.24 

2.60+0.19 

- 


19014228-0027388 

5.16+0.18 

- 

5.86+0.20 

6.51+0.25 

2.35+0.14 

- 


19014498-0027496 

5.12+0.11 

- 

5.82+0.16 

6.37+0.29 

2.33+0.22 

- 


19014525-0023580 

5.17+0.11 

4.18+0.12 

5.84+0.17 

6.46+0.22 

2.44+0.28 

0.69+0.06 


19014769-0025108 

5.17+0.13 

4.15+ 0.17 

5.81+0.20 

6.39+0.24 

2.54+0.11 

0.67+0.10 


19015261-0025318 

5.17+0.10 

- 

5.91+0.16 

6.46+0.25 

2.53+0.19 

- 


19015978-0028183 

4.74+0.10 

3.65+ 0.22 

5.15+0.13 

5.81+0.18 

2.32+0.26 

0.39+0.27 



Table 6. Results, errors, and estimated systematic uncertainties using different evolutionary models with [Fe/H]~0.2. 


Model 

age 

(Gyr) 

(m - M)o 
(mag) 

E(B - V) 
(mag) 

do 

(kpc) 

Rcc 

(kpc) 

z 

(pc) 

Mto 

(Mo) 

PARSEC 

0.98+0.1 

12.72+0.1 

0.83+0.02 

3.50 

5.45 

-152.18 

2.2 

BASTI 

0.75+0.1 

12.85+0.1 

0.86+0.02 

3.72 

5.33 

-161.57 

2.3 

BASTE 

0.65+0.1 

12.75+0.1 

0.89+0.02 

3.55 

5.42 

-154.30 

2.3 

Dartmouth 

0.90+0.15 

12.72+0.1 

0.84+0.04 

3.50 

5.45 

-152.18 

2.2 

Dartmouth" 

0.80+0.15 

12.70+0.1 

0.83+0.04 

3.47 

5.47 

-150.79 

2.2 

Victoria-Regina 

0.90+0.15 

12.75+0.1 

0.85+0.04 

3.50 

5.42 

-154.30 

2.3 


"a-enhanced models: [a/Fe]—1-0.4 dex in the case of BASTI (we chose the models with [Fe/H]=0.16 dex, the closest metallicity 
to the spectroscopic derivation) and [a/Fe]=-1-0.2 dex ([Fe/H]=0.23 dex) in the case of Dartmouth. 

Table 7. Clusters’ parameters 


Name 

RA dec 

J2000.0 

1 

deg. 

b 

deg. 

E{B - V) 

Age 

(Gyr) 

RGc(a) 

(kpc) 

Ref. 

Be 81 

19:01:36 

-00:31:00 

34.505 

-2.07 

0.85 

0.98±0.1 

5.45 

present work 

NGC 6705 

18:51:05 

-06:16:12 

27.31 

-2.78 

0.43 

0.30±0.05 

6.3 

idrI, Cantat-Gaudin et al. (2014) 

Tr20 

12:39:32 

-60:37:36 

301.48 

+2.22 

0.33 

1.50±0.15 

6.88 

idrI, Donati et al. (2014b) 

NGC 4815 

12:57:59 

-64:57:36 

303.63 

-2.10 

0.72 

0.57±0.07 

6.9 

idrI, Friel et al. (2014) 


(a) computed with Rq= 8 kpc. 
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Table 8. Stellar Nucleosynthesis 


Element 

Main Production Site 

MECHANISM 

yield(SN Ia/SN II) 

IdQ 

Massive Stars 

Helium burning 

8% 

“Mg 

Massive Stars 

C, Ne burnings 

10% 

“Si 

Massive Stars 

explosive and non-explosive 0 burning 

60% 

®Ca 

Massive Stars 

explosive and non-explosive 0 burning 

67% 

«Sc 

Massive Stars 

C, Ne burnings, a and v-wind (neutrino-powered wind) 

49% 

“Ti 

Massive Stars and SNIa 

explosive Si burning and SNIa with He detonation 

63% 

5‘V 

Massive Stars and SNIa 

explosive Si and 0 burnings, SNIa with He detonation, and a and v 

88% 


Massive Stars and SNIa 

explosive Si burning, SNIa with He detonation, and a 

84% 

“Fe 

Massive Stars and SNIa 

explosive Si burning and SNIa 

88% 

“Ni 

Massive Stars (and SNIa) 

or (o'-rich freeze-out from nuclear statistical equilibrium) and SNIa 

75% 

SOy 

Massive Stars 

He-burning s-process, and v-wind 

- 


Massive Stars, compact binary merger 

v-wind 

- 
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